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Abstract 

A soluble model of weakly coupled "molecular" and "nuclear" Hamiltoni- 
ans is studied in order to exhibit explicitly the mechanism leading to the 
enhancement of fusion probability in case of a narrow near-threshold nuclear 
resonance. We, further, consider molecular cells of this type being arranged 
in lattice structures. It is shown that if the real part of the narrow nuclear 
resonance lies within the molecular band generated by the intercellular in- 
teraction, an enhancement, proportional to the inverse width of the nuclear 
resonance, is to be expected. 

PACS numbers: 03.65.-w, 05.50.-hq, 24.30.-v, 25.60.Pj 

1. INTRODUCTION 

Molecules are usually treated as purely Coulombic systems, while the strong interaction 
between their nuclear constituents is assumed to play a negligible role. However, at least 
in principle, any Coulombic molecular level lying above the lower threshold of the nuclear 
subsystem, is embedded in the continuous spectrum of the nuclear sub-Hamiltonian. The 
coupling between the molecular and nuclear channels, hence, turns this level into a resonance 
(see, e. g., Refs. [TJ-Q and references cited therein). Of course, due to the wide Coulombic 
barrier between the nuclei and the short-range character of the nuclear interaction, this 
coupling, and thus the width of the resonance, which determines the fusion probability of 
the nuclear constituents of the molecule, is in general extremely small. 

However, as pointed out in [^,[7| , the situation is rather different if the nuclear subsystem 
of a molecule has a sufficiently narrow near-threshold. Examples of such nuclear systems 
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may be read off from the data presented in Among them are even customary systems 
hke pp^^O and p^^O |^,10|, i.e., the nuclear constituents of the water molecule H2O or 



the hydroxyl ion 0H~ with O being the isotope ^^O. For LiD and H2O the influence of 
near-threshold nuclear resonances on the molecular properties has been studied in |1TI|-|T3| by 
estimating the overlap integrals between the corresponding molecular and nuclear wave func- 
tions. The best known example of such phenomena is the muon catalyzed fusion of deuteron 
and triton in the dtfi molecule, where the near-threshold nuclear resonance ^He(3/2"'") plays 
a decisive role [|1^]. 

Being motivated by the above special cases, we deal in this paper with a rather general 
model Hamiltonian related to the ones considered by Friedrichs in [|15|. This Hamiltonian 
consists of a "nuclear" part, a "molecular" part with eigenvalues embedded in the continu- 
ous spectrum of the nuclear part, and a weak coupling term which turns these unperturbed 
eigenvalues into "molecular" resonances. Since the model is explicitly solvable, the mecha- 
nism of formation of the resonances becomes clearly visible. 

The following property pointed out in appears, in particular, as a general feature: if 
the "nuclear" channel itself has a narrow resonance with a position close to the "molecular" 
energy, then the width (the imaginary part) of the resulting "molecular" resonance is found 
to be inversely proportional to the "nuclear" width. In other words, a large increase of the 
decay rate of the "molecular" state, i. e. of the fusion probability, is observed in this case. 
Such a coincidence of nuclear and molecular energies is, of course, a rather rare phenomenon 
in nature. 

A further goal of the present work is to show that the decay rate may be considerably 
enhanced when arranging molecular clusters of this type within a crystalline structure. The 
reason is that in such a configuration the original discrete molecular energy turns into a band, 
i.e., into a whole interval of the continuous spectrumQ. That is, even if the position of the 
"nuclear" resonance differs from the original "molecular" level, it can get into this band. This 
allows for a fine tuning by exciting the crystalline structure to energies as close as possible 
to the energy of the "nuclear" resonance. We show that the lattice states, which correspond 
to such an initial choice of their quasimomentum distribution, decay exponentially with a 
rate which is again inversely proportional to the width of the "nuclear" resonance. 

The paper is organized as follows. 

In Sec. we introduce the explicitly soluble model designed to demonstrate the interplay 
of the molecular and nuclear resonance widths. It is also shown that in a wide time interval 
the decay of the "molecular" state is indeed of the standard exponential character WE 



This transition will take place primarily into the open nuclear channels and its rate is 
determined by the inverse width of the nuclear resonance. In Sec. ^ we consider the case 
where molecular Hamiltonians of the type considered in Sec. |^ are arranged in form of an 
infinite one-dimensional lattice. Sec. § is devoted to the generalization to mult i- dimensional 
lattices. In both these sections, the time evolution of originally pure molecular states, 
extended to a spectral band within a lattice, is studied. It is shown that if the real part of 



Note that in the models under consideration, the spectral band generated by the "molecular' 
level is shifted finally, after switching on the coupling between the "nuclear" and "molecular' 
channels, into the unphysical sheet. 
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the "nuclear" resonance lies within such a spectral band, then there exist molecular states 
which decay exponentially, with a rate inversely proportional to the "nuclear" width. 



2. TWO-CHANNEL MOLECULAR RESONANCE MODEL 

A. Description of the model Hamiltonian 

Let us consider a two-channel Hilbert space Ti = Tii ® 7^2 consisting of a "nuclear" 
Hilbert space Tii (channel 1) and a one- dimensional "molecular" space 7^2 = C (channel 2). 

Ml 
U2 

U2 being a complex number. The inner product {u,v)h = {ui,vi) + U2V2 in 7i is naturally 
defined via the inner products {ui,vi) in Hi and U2V2 in 7^2- 
As a Hamiltonian in Ti we consider the 2 x 2 operator matrix 



The elements of 7i are represented as vectors u 



where Ui G Hi and U2 &H2, with 



A 



hi b 
{■,b) A2 



(2.1) 



where hi is the (selfadjoint) "nuclear Hamiltonian" in Hi, and A2 € M a trial "molecular" 
energy. A vector b eHi provides the coupling between the channels. It should be mentioned 
that the Hamiltonian (|2.1| ) resembles one of the well known Friedrichs models [1^. 

If there is no coupling between the channels, i. e. for 6 = 0, the spectrum of A consists of 
the spectrum of hi and the additional discrete eigenvalue A2. We assume that the continuous 
spectrum (Tc(/ii) of the Hamiltonian hi is not empty and that the eigenvalue A2 is embedded 
in (Tc(/ii). It is also assumed that A2 is not a threshold point of crc(^i), and that this spectrum 
is absolutely continuous in a sufficiently wide neighborhood of A2. 

A nontrivial coupling between the channels will, in general, shift the eigenvalue A2 into an 
unphysical sheet of the energy plane. The resulting perturbed energy appears as a resonance, 
i. e., as a pole of the analytic continuation of the resolvent r{z) = {A — z)~^ taken between 
suitable states (see, e. g., 0). The explicit representation for the resolvent r{z) is easily seen 
to be 



r{z) 



ri{z) 



ri{z)b{ ■,b)ri{z) ri{z)b \ 



\ 



M2{z) 
{ ■ ,b)ri{z) 

M2{z) 



M2{z) 

1 

M2{z) J 



(2.2) 



Here ri{z) stands for the resolvent ri{z) = {hi — z)^^ of the "nuclear" Hamiltonian hi, while 
the transfer function M2{z) reads M2{z) = X2 — z — j3{z) with l3{z) = {ri{z)b, b). Evidently 
the poles of r{z) on the physical sheet are either due to zeros of the function M2{z) or due to 
poles of the resolvent ri{z). The latter correspond to the discrete spectrum of the operator 
hi which may determine part of the point spectrum of A. This is true, in particular, for the 
multiple eigenvalues of hi. In any case it is obvious that the perturbation of the eigenvalue 
A2 only corresponds to solutions of the equation M2{z) = 0, i.e., of 



z = X2-(3{z)- 



(2.3) 
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This equation has no roots z with Im^; 7^ on the physical sheet. For, being eigenvalues 



of the selfadjoint operator A, they have, of course, to be real. Thus, Eq. (2^) may have 
solutions only on the real axis and in the unphysical sheet (s) of the Riemann surface of the 
resolvent rxi^z). 

We start with a brief discussion of the case where the "nuclear" channel Hamiltonian 
h\ generates no resonances close to A2 in a domain T> of the unphysical sheet which ajoins 
the physical sheet from below the cut. This assumption implies that for a wide set of unit 
vectors h = b/\\b\\ the quadratic form P{z) = \\b\\'^{ri{z)b, b) can be analytically continued in 

I, Eq. ( |2.3| ) is uniquely solvable ||n 
see, e.g., 



T). Moreover, under certain smallness conditions for 
in V providing in first order perturbation theory 



Z2 = A2-(ri(A2 + iO)6,6)+o(||6|n. 



(2.4) 



The real and imaginary parts of the resonance Z2 = E^^ — i— thus, are given by 



r 



(2) 



E 



(2) 
R 
.(2) 
R 



A2 -Re(ri(A2 + iO)6, b) 
2Im(ri(A2 + iO)6, b) + o 



(2.5) 



B. Perturbation of the "molecular" resonance by a nearby "nuclear" resonance 

Our main interest concerns the opposite case of a "nuclear" resonance zi = E^j^^ — 

r^-* > 0, with a real part E^j^^ close to A2. For the sake of simplicity we assume the 
corresponding pole of ri{z) to be of first order. Let the element 6 G TYi be such that the 
function (3{z) admits an analytic continuation into a domain V which contains both points 
A2 and Zi. This domain, moreover, is assumed to belong to the unphysical sheet which 
adjoins the physical sheet along the upper rim of the cut. In T> the function /3(z), thus, can 
be written as 

Piz) = ^— + r^{z) (2.6) 
Zi — z 

with (3'^^^{z) being a holomorphic function. For a fixed "structure function" b = b/\\b\\ we 
have \a\ = Ca||6||^ with a constant Ca determined by the residue of ri{z) at z = z\. Note 
that this residue is usually expressed in terms of resonance (Gamow) functions (see for 
example pO|)- ■'■'^ fact, we assume that the resonance corresponds to an "almost eigenstate" 



of h\. That is, in principle a limiting procedure F^'' — > is possible so that the resonance 
turns into a usual eigenvalue with an eigenvector G Tii. More precisely, we assume 

= Cf + o(l) as F«^0 (2.7) 

with = ^ 0. This can be achieved, e. g., if the Hamiltonian hx itself has a 

matrix representation of the form ( p.l| ) and the resonance Z\ is generated by a separated one- 



dimensional channel. In such a case, according to (|2.2|) and ( p.4|) , we would have Cf^ = 1 
(for details see Ref. 0, Sec. II). 
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Let 



Re a > and Im a -C Re a 



(2.8) 



and, for z eV, 



lmP'''^{z)\>cr,\ 



and \P''^iz)\ <Cv\\bf . 



with constants Cx) > and > 0. Furthermore, the couphng between the channels in the 
Hamiltonian (U-l\) is assumed to be so weak that 



\r^{z)\<G 



x)||y|P ^ r^-* while \a\ 



Ca\ 



< r 



^(1) 

R 



(2.9) 



It can be expected that these conditions are fulfilled in specific molecular systems even under 
the supposition that the "nuclear" width F^'' itself is very small. 

After inserting (^]6|) for P{z), Eq. ( p. 3D turns into the "quadratic" equation 

(A2 -z){z,-z)-a+ (zi - z)(3''^{z) = 

which can be "solved", i.e., can be rewritten in form of two equations 



± 



+ a. 



(2.10) 



Under conditions (|2.9| ) the existence of solutions of ( 2.1C ), and thus of Eq. (|2.3| ), is guar- 
anteed, analogously to the proof in |]l9l, by Banach's Fixed Point Theorem. Each of the 
equations (p.lO|) has only one solution in the domain V. In case of the sign "— " we denote 
the root of (|2.10|) by 2;nuci, in case of the sign "+" by Zmoi- 
The inequalities (p.9[) imply 



X9.-E 



(1) 



+ 



F 



(1)" 



< 



4CJ|fc| 



F 



(1) 



< 1. 



(2.11) 



For z G the value of 



4a 



[X,-z,~l3^-^{z)Y 



(2.12) 



is very small, \b{z)\ ^ 1. Thus, to separate the main terms of the solutions of Eqs. ( [2.1(J ), 
one can apply the asymptotic relation y/1 + e = 1 + £:/2 + 0{e'^). As a result we find 



2a 



2 2 V (A2-^i-/3-K^))2 



In other words, the roots z^^ci and z^oi of ( |2.10| ) are essentially given by 

a a 



'2'nucl — Zi 



X2-Z1- /?-s(zi) 

Zmol = X2-r^{X2 + iO) + 



Zl 



X2 - Zl 
a 



X2-Z1- /?'^<=s(A2 + iO) 



A2 + 



Zl 



(2.13) 

(2.14) 
(2.15) 
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^From the second condition (|2.9| ) follows 



a 



^ r^-* . Consequently, this term provides 



\2- Zi 

in 2;nuci a very small perturbation of the initial "nuclear" resonance Zi. As compared to F^'' 
it represents also in z^ox a very weak perturbation of the "molecular" energy A2. However, 
as compared to the result (p.4|), valid in case of a missing nearby "nuclear" resonance, it can 



be rather large. In particular, if the "molecular" energy A2 coincides with the real part e'^^ 



of the "nuclear" resonance Zi, then Zmoi = e'^r^ — i ^ with 



(m) 



4'"'^A,-2^ aud rr^45i^, (2,16) 



^ R ^ R 

The width of the "molecular" resonance z^oi in the presence of a nearby "nuclear" resonance 
Z\, thus, turns out to be inversely proportional to the "nuclear" width T^^\ 

Let us contrast the results ( |2.4| ) and ( |2.16| ) in some more detail. Since such a comparison 
is necessarily of a qualitative character, we simulate the situation of a missing nearby nu- 
clear resonance simply by dropping the pole term in the representation ( |2.6| ) of l3{z). After 
this removal we get f3{z) = l3'^'^^{z) and for Im^ < the eigenvalue A2 generates the reso- 
nance (|2.4| ) having the width F^"* ^ 2 Im/3'''^^(A2 + iO). The latter satisfies the inequalities 

cv\\b\\^ < FgV2 < Cv\\bf. Substituting |Rea| ~ Ca\\bf ~ — Fg^ in (gJ§ we find the 



following approximate estimate of F^-* relative to Fg'*: 



Fi™) ~ Fg) . (2.17) 

^ R 

The second inequality ( p.9| ), chosen as a condition for reflects the fact that the "usual" 
molecular width Fg'' is much smaller than the width of a usual "nuclear" resonance Fg'*, 

CargU<c,,(Fg))'. (2.18) 

This can practically always be assumed for concrete molecules. 

Under condition ( ^TTD the value of Ca = differs from zero, > C > 0, as 

F^"* —>■ 0. Therefore the estimates ( p. 16 ) and ( p. 17 ) imply that in the presence of a narrow 



Ca/Cv 



[F^-* <^ Ca/cxi) "nuclear" resonance close to A2 the "molecular" width F^"'* is much larger 

such a 

Ca/cv 



(2) 

than the "molecular" width V)^ observed in absence of such a resonance. In fact, according 



to ( |2.17 ) this ratio is determined by the large quotient 



r(i) 

^ R 



Finally we note that if the conditions ( |2.9|) , and hence the condition (|2.18 ) are not 



fulfilled, i.e., if the coupling between the channels in the Hamiltonian ( p.l| ) is not small 
compared with the "nuclear" width, then it follows from ( p.lOj ) that the molecular width 
Fg^'' achieves itself an order of . We do not discuss this case since such a situation can 
hardly be assumed to exist. 
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FIG. 1. A scheme showing the deformation of the integration path 7. The part 7"'' of the 
resulting contour belongs to the unphysical sheet, part 7" to the physical sheet. The crosses 
"x" denote the discrete eigenvalues of the Hamiltonian A while the solid line corresponds to the 
continuous spectrum. 

C. Exponential decay of the "molecular" state 

Let us suppose that an initial state of the system described by the Hamiltonian A corre- 
sponds exactly to the pure "molecular" wave function cp = ( ? ) . Then, the time evolution 



^ 1 ^ 

of the system is described by the solution ip{t) of the Cauchy problem 

^^=M, ^\t=o = ^- (2-19) 
The probability of finding the system at the time t still in the molecular state <f is given by 

Praol{t) = \{m,f)\'- 

The remainder 1 — Pmo\{t), hence, determines the probability for the state ip to decay into 
open channels of the continuous spectrum of the "nuclear" sub-Hamiltonian hi. 

To estimate the probability Pmoi(^), we use the standard integral representation of a 
function of an operator via its resolvent. In the case considered this means 

exp{-iAt} = j dze~''''\A- z)-^. (2.20) 
7 

The integration in (|2.20|) is performed in the physical sheet along a contour 7 going counter- 
clockwise around the spectrum of the matrix A. Recall that, due to the selfadjointness of the 
operator A, this spectrum is real. Taking into account the representations ( p.2|) and (p. 20 ) 
one finds 



This leads to the following important result. Under the conditions of Subsection |2 B| the 
behavior of the integral (|2.21|) for t > is described by the formula 



7 



{i'(t),(p) = exp{-i2;moit} 



+ exp{-i2;nuci0 



A2 - ^1 - /9'-'=g(A2 + iO))" 



+ o(£^(A2 + iO)) 



(2.22) 



(a2-;2i-/?-s(^i))' 



where the value of e{z) is given by Eq. ( p.l2| ). The background term 6{t) = 0(||6p) is small, 



\£{t)\ <^ 1, for all t > 0. In particular, if E'^-' = A2 we have 



{i){t),ip) ^ exp{-i2;„ioit} 



4a 



^(1) 

R 



CXp{-iZnuclt} 



4a 



.(1) 

R 



+ 



(2.23) 



The proof of the asymptotic relation ( |2.22| ) is carried out by estimating the contribution 
of the resonance poles z^o\ and z^^^ci to the integral (|2.21|) . This is done by deforming parts 
of the contour 7 situated in a neighborhood of the molecular energy A2 (see Figure |I]). A 
part 7"*" of 7, situated initially on the upper rim of the cut, is shifted into the neighboring 
unphysical sheet. Having done such a deformation one finds explicitly the residues of the 
integrand in ( |2.21| ) aX z = z^ox and z = ^nuci- An analogous deformation of a part 7" of 
7, situated initially on the lower rim, is performed in a domain Im z < of the physical 
sheet. It is assumed that, though the parts 7"*" and 7" belong to different energy sheets, 
their positions on these sheets coincide. It is also assumed that for any z G 7^ the estimate 
1/3^(2;) I ^ IA2 — z\ holds. Thus, the integration in ( [^.21| ) around the continuous spectrum 
of A, except the residues at z = z^^oi and z = ^nuci, gives 



-— yd^exp( ,zt} [^^_^_^+(^)][^^_^_^-(,)]- 



(2.24) 



Here we have used the specific notation (3~^{z) for the values of the function (3{z) at points 
z belonging to the curve 7"*" (i. e., lying in the unphysical sheet), and P~{z) for the values of 
P{z) at the same points of the curve 7" (i.e., lying in the physical sheet). Both P~{z) and 
P~^{z) are of the order of 0(||6|p), and |/9^(z)| -C IA2 — z\, while the exponential exp{—izt) 
at Imz < is decreasing for t > 0. The value of the function ( p.24| ), thus, is always 
small, having an order of 0(||6|p), and is even decreasing (in general nonexponentially) 
with increasing t. We include the contribution of this function in the background term 
£{t). The summand e{t) also includes a contribution to (|2.21 ) from the residues at the 
discrete eigenvalues of A. Apart from factors oscillating when t changes, the value of this 
contribution remains practically the same for all t > 0. 

The formulae ( p.22| ) and ( p.23| ) show explicitly that in a large time interval < t < T, 

T ~ ~(mj liiKiax |e(t)| , the decay of a "molecular" state if in the presence of a narrow 

^R 
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"nuclear" resonance is indeed of exponential character. The rate of this decay is determined 
mainly by the width F^"* of the "molecular" resonance ^moh by the ratio | Rea|/r^'', 

P^oiW = exp{-T^^h} = exp <| -^^^ t \ . (2.25) 



'MOLECULAR" RESONANCES IN A ONE-DIMENSIONAL LATTICE 



Let us assume that the "molecules" described by the Hamiltonian are arranged in 
form of an infinite one-dimensional (linear) crystalline structure. To describe such a crystal 
we introduce the lattice Hilbert space 

+00 

g= ® n^'^ 

i=—oo 

representing an orthogonal sum of the Hilbert spaces associated with the individual cells 

=1-Lf ®1-Lf . (3.1) 

Here the subspaces T-^i = Hi and 7^2^ = 7^2 = C are exactly the same ones as in Sec. ^ and, 
thus, Ti^*-* = 7i. The elements of the total Hilbert space Q are represented by the sequences 

u = {...,u^ '^\u^ '^\u^^\u^^\u'^'^\ . . .) with components -u^*^ = L ) where u{ E Hi 
and U2 E H2 = C The inner product in H is defined by {u,v)u = Yl {u^^\v^^^)'h(^)- The 

i=— 00 

+00 +00 /j-, 

subspaces ^1 = (B Hi and Q2 = © H2 , with Q = ^1 © ^2, are called pure nuclear and 

i=— 00 i=—oo 

pure molecular channels, respectively. 

In the present section we deal with the Hamiltonian H acting in H according to 

(//m)» = Wu^'-'^ + Au^''^ + Wu^'^'^ (3.2) 

where only the interaction between neighboring cells is taken into account and the interaction 
operator W is chosen in the simplest form 

:) (3-3) 

with w being a positive number. Such a choice of the interaction corresponds to the natural 
assumption that the cells interact between each other via the molecular states, while the 
direct interaction between nuclear constituents belonging to different cells is negligible. We 
assume that the closed interval [A2 — 2w, A2 + 2w] is totally embedded in the continuous 
spectrum o"c(/io) of ^0 smd, moreover, that no thresholds of ac{ho) belong to this interval. For 
the sake of simplicity we also assume that the interval belongs to the domain V introduced 
in Sec. ^ and that for any E [X2 — 2w, A2 + 2w] 
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Im(ro(/i±iO)6,6) ^ 0. (3.4) 

Obviously, the Hamiltonian ( p.2| ) represents a special case of the infinite Jacobi operator 
matrix (regarding the properties of some infinite scalar Jacobi matrices see, e.g., |]2l|j2^ 



and Refs. cited therein). It is a selfadjoint operator on the domain T>{H) = © V^''' with 

i=— oo 

p« = V{hi)®C. 

The resolvent R{z) = {H — z)^^ of H possesses a natural block structure, R{z) = 
{R{j, k] z)}, j, k = 0, ±1, ±2, ±oo. The blocks R{j, k; z), providing the mapping 
'H^^\ satisfy the equations 

WR{j - 1, k; z) + {A- z)R{j, k; z) + WR^j + 1, k; z) = Sjkl (3.5) 

where 6jk stands for the Kronecker delta and / for the identity operator in the Hilbert 
space H of cells. Hereafter we assume Im^ 7^ so that the value of z automatically 
belongs to the resolvent set of the operator H. The blocks R{j, k; z) themselves possess 
a 2 X 2 matrix structure, R{j,k]z) = {Rmn{j, k; z)}, m,n = 1,2, corresponding to the 
decomposition H = Hi (B'H2- 
The Fourier transform 

^ +00 

{Fu){p) = — == V u^^^ exp(ipj) (3.6) 
V Ztt — 

j=-co 

in Q reduces Eq. ( p.5| ) to 

{A - z)R{p, p'; z) + 2 cospWR{p, p'; z) = S{p - p') I (3.7) 

where the quasi-momentum p runs through the interval [— vr, vr] and the function R[p,p'] z) 
represents the kernel of the resolvent R{z) in this representation. From ( p.7|) follows imme- 
diately 

Rip,p';z)=Gip;z)5{p-p') (3.8) 

where 



G(p; z) 



( / N , rx{z)h{- ,h)rx{z) ri{z)b \ 

I ri{z) H — -— » 

M2{p;z) M2{p;z) 

{■,b)niz) 1 



(3.9) 



V M2{p;z) M2{p;z) / 



This corresponds to the representation ( |2.2| ) of the resolvent of the cell Hamiltonian A, the 
only difference being that the transfer function M2{z) is replaced by the expression 

M2{p;z) = X2- z + 2w cos p- (3 (z) . (3.10) 
The factorization ( p.8| ) implies 

R{j,k;z) = ^ r dpe-'P^^''^G{p;z) 
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and with the representation ( p.9|) we, thus, obtain 

JD(■t..^^_^ Sjkri{z) + ri{z)bR22{j,k;z){- ,b)ri{z) -ri{z)b R22{j, k; z) \ 

' ^ V ~R22{j,k;z){-Sri{z) R22{j,k;z) J ^^'^'^ 

where 

R22{j,k;z) = — dp- — -—. (3.12) 

zn X2 — z + 2w cosp — p{z) 

Introducing the new variable ( = e~^^, this integral is reduced to 

X r 

R22U,k;z) = — d>d( 



2m J ^ wC + M2{z)C + w 
7 

Here, 7 stands for the unit circle centered at the origin, the integration over 7 being per- 
formed in the counterclockwise sense. Further, applying the Residue Theorem and taking 
into account the sign of the difference j — k one finds 

^/[M2{z)-2w][M2{z) + 2w] - M2{z) 



R22U, k; z) = ^^^^ , . (3.13) 

^[M2{z)-2w][M2{z) + 2w] 

It is assumed here that the branch a/(^ — 2w)(^ + 2w) of the function ^(^ — 2w)(^ + 
\ 1/2 

2w)] is defined in the plane of the complex parameter ^, cut along the interval 



[— 2w, 2w], and that Im^ > implies Im \/ — 2w)(^ + 2w) > 0, while Im^ < implies 
Im ^(^-2w)(^ + 2w) < 0. 

iFiom. Eqs. (|3.11| ) and (|3.13|) it follows that all the nontrivial singularities of the resolvent 
R{z), differing from those of the cell "nuclear" channel resolvent ri{z), are determined by 
the properties of the function 

D{z) = ^[M2{z) -2w][M2{z) + 2w]. 

First, we note that if = 0, and thus M2{z) = X2 — z, then the "molecular" and "nuclear" 
channels in the Hamiltonian H decouple. In this case the eigenvalue A2 generates for H 
an additional branch of the absolutely continuous spectrum which occupies the interval 
[A2 — 2w, A + 2w]. Second, even if 7^ then the function D{z) cannot have roots z 
with Im2; 7^ in the physical sheet. Otherwise such roots would generate for H a complex 
spectrum. But this is impossible because of the selfadjointness of H . Also, under the 
condition (|3.4| ) this function cannot have real roots within the interval [A2 — 2w, A2 + 2w] 
since for A2 — 2w < < A2 + 2w the imaginary part 

Im[M2(/i ± iO) - 2w] = Im[M2(/i ± iO) + 2w] = -\\bf Im(ro(/i ± iO)6, 6) 

is nonzero. Thus, in a close neighborhood of the interval [A2 — 2w, A2 + 2w] the equation 
D{z) = may only have roots in the unphysical sheet. In fact, assuming the conditions ( p.9|) 
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X XX 



A2 — 2w A2 A2 + 2w 
1 ■ \ 




FIG. 2. A scheme showing the position of the resonance bands generated in the unphysical 
sheet by the "molecular" eigenvalue A2 and the "nuclear" resonance zi . These bands are generated 
respectively by the points z^oiip) a-nd z^^^iip) with the quasimomentum p running through the 
interval [— tt, vr]. 



and repeating literally the considerations which led to ( p.l3|) , one can rewrite this relation 
in form of the four equations, 



A2 — 2w + Zl — P 



reg/ 



± 



A2 - 2w - Zl - p'^^iz) 



X2 + 2w + Zl - P'^^jz) 
2 

X2 + 2W-Z1- P'^^^iz) 



1 + 



2a 



(A2 -2w-zi- p'^^iz)y 



(3.14) 



(3.15) 



1 + 



2a 



(A2 + 2w- zi-p'' 



where e± = 4a/(A2 ± 2w — Zi— P'^'^^{z)). In that part of the domain V which belongs to the 
unphysical sheet, equation D{z) = has four solutions being given essentially by 



,(±) 

"nucl 
^mol 



a 



a 



Zl - 



A2 ± 2w - Zl - /3''^s(zi) 
A2 ± 2w - /3''°^(A2 ±2w + iO) + 
a 



Zl 



A2 ± 2w — Zl 
a 



A2 ± 2w - ^1 - /3'-^s(A2 ±2w + iO) 



(3.16) 
(3.17) 



A2 ± 2w + 



A2 ± 2w — Zl 



Obviously, each of the roots z^^^-y and represents an additional square-root branching 
point of the Riemann surface of the functions R22{j, k; z). Consequently, these roots are 
also the branching points of the Riemann surface of the total resolvent R{z). Thus, one 
has to introduce the "resonance" cuts in the unphysical sheet considered. The cuts can be 
made, say, between z^~^\ to z^^l-^ and between 2;^] and z^|. Evidently, these cuts are to 
be interpreted as the resonance spectral bands generated by the initial "molecular" level A2 
and the "nuclear" resonance zi (see Figure 0). 

Consider now the time evolution of the system described by the Hamiltonian H starting 
from a pure molecular state 93 = yji © ip2, \\^Pn\\ ^ Gn, n = 1,2, with ipi = and \\{p\\ = 
1 1^32 1 1 = 1- The probability to find the system at a time t > in the molecular channel is 
given by 
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i-fft,_l|2 



(3.18) 



where P2 is the orthogonal projection in Q on the pure molecular subspace Q2- 

As in (|2.2CI| ) we represent the time evolution operator exp(— lift) in terms of the resolvent 

R{z) = {H-z)-\ 



exp {-iHt} = --^T (f) dze-''\H - z)'^ 

ZTTl 



(3.19) 



where the integration is performed along a counterclockwise contour 7 in the physical sheet 
encircling the spectrum of the Hamiltonian H. Recall that this spectrum is real since H is 
a selfadjoint operator. 

According to Eqs. ( p.8|) and ( p.lO|) the operator P2R{z)\g^ acts in quasi-momentum 
representation as the multiplication operator 



G22ip,z) = [M2{p; z)]-' 



(3.20) 



and, thus, 



P2R{Z)^){P) 



1 



Mp) , P e [-7r,7r] . 



M2(p; z) 

Here (p2{p) stands for the values of the Fourier transform ( p.6|) of the vector (p2 

, (-2) (-1) (0) (1) (2) s , . , 

{■■■■,^2 \'^2 ,V52 ^V^2 5 V52 ^ which means 



P2e-'^V)(P) = -^ ^2{p) J{p,t) 



with 



J{p,t) = (p dz 



exp(— izt) 



X,{p) - z - /3iz) 



(3.21) 



(3.22) 



This expression has exactly the same form as the integral (|2.21 ) . The only difference consists 
in the replacement of A2 by the sum 



X2{p) = A2 + 2w cosp . 



(3.23) 



Thus, to estimate the function J(p, t) one can immediately use the relation (|2.22|) in order 
to find 



J{p,t) = exp{-iZraoi{p)t} 



+ exp{-iznuci(p)0 



A2(p)-^i-/?-g(A2(p) + iO) 



+ 0{6\p,\2ip) + iO) 



A2(p) -zi- P'^^izi) 



0{s\p,z,] 



+ e{p,t), (3.24) 



where 
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eip, z) = ^ . (3.25) 

[A2(p)-^i-/3-g(^)]2 



The function e{p,t) = 0(||6|p) is always small, t)| ^ 1. In accordance with Eqs. (|2.14|) 
and (|2.15|) we, hence, obtain for the positions of the poles 



^nucl(p) = z,- —— , (3.26) 

A2 + 2w cosp — Zi 

2;moi(p) = A2 + 2wcosp + - — . (3.27) 

A2 + 2w C0SJ9 — Zi 

The resonance bands representing the ranges of the functions Znuci(p) and ^moi(p), with p 
running through the interval [— 7r,-7r], are schematically depicted in Figure |[ 
The asymptotics ( p.24| ) implies 

Praoi{^,t)= r dp\Jip,t)\'\Mp)\'= r dpexp{-T^^\p)t}\Mp)\' + m (3.28) 

J —TV J —TV 



where 



T^^\p) = -2lmz^M = -2Im . " • (3.29) 

A2 + 2wcosp — 2:1 



The background term e{t) in ( |3.28D is small for any t > 0, e{t) = 0(||6|| ) and \e{t)\ <^ 1. 



Further, let us assume that the real part E^j^^ of the "nuclear" resonance zi belongs to 
the interval [A2 — 2w, A2 + 2w], that is l-^jj^^ — A2I < 2w. Then, one can always prepare an 
initial "molecular" state ip which decays via the "nuclear" channel with a rate as close as 
possible to the desired maximal value. Under the assumption ( p. 81 ), this maximum is given 
by 

max FV, = A—rrr- 

^ R 

(cf. formula ( p.l6|) ). The correspondingly prepared "molecular" state (p has an almost 
monochromatic component 'p>2{p) being localized in a close neighborhood of the quasi- 
momenta 

I ^R ~ ^2 

» = ± arccos — ^ . 

^ 2w 

For example, if the function ^2{p) is nonzero only for quasimomenta p restricted by 



Eg) - A2 

cos p — 



2w 



< 5 — ^ with some small 5 > 0, then the width Fp given by the rela- 

- Aw R s J 

___ . . . . 1 4 Re a , 4 Re a 

tion ( p.29|) varies m an mterval lymg approximately between — — — and — jjy. 

R R 

4. "MOLECULAR" RESONANCES IN A MULTI-DIMENSIONAL LATTICE 

In this section we consider the case where the "molecules" described by the Hamil- 
tonians (|2.1| ) form an infinite A^-dimensional crystalline structure. To label the cells 
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of the respective lattice we use the multi-index i G Z^, i.e., i = {11,12, ...^in) with 
ik = ...,—2,-1,0,1,2,..., k = 1,2,...,N. The Hilbert space of the system considered 
is in this case Q = (B 'H^^\ where the individual cell spaces are given by (|3.1j), with 

Hi^ = Hi and 7^2'' = H2 = C being the spaces introduced in Sec. |^. For the components 
M^*^ G Ti*-*-* of the elements u of the total Hilbert space Q we again use the column repre- 

sentation -u*^*) = ( L ) with u^^ G Hi and & H2 = C The inner product in H is 

defined as {u,v)'h = ^ f (*^)-^{,). The subspaces ^1 = © Hi^ and Q2 = © '^2'') 

with Qi (B G2 = G, represent pure nuclear and pure molecular channels, respectively. 
The Hamiltonian is defined in Q by the expression 

(i^M)» = AM»+ J2 Wii,j)u^^^ (4.1) 



where the interaction matrices 



Wll[i,J) Wi2{l,J, 
W2lii,j) W22ii,j] 



consist of the block components Wmnihj) providing the mappings Hn — > Hm, m,n = 1,2. 
These components describe the direct interaction between the m-th and n-th channels 
of the different cells i and j, respectively. The matrices W{i,j) are assumed to be 
bounded operators in H which depend only on the difference i — j = {ii — Ji,^2 — 
j2,...,iN — Jn), i-e. W{i,j) = W{i — j) and, thus, the same holds for the block com- 
ponents, Wmnihj) = Wmn{i ~ j)- Moreover, the series of W{j) is assumed to be convergent 
with respect to the operator norm topology, i. e. 

J2 ||W^(j)ll<+oo (4.2) 

and the property 

Wij-t) = [Wit-j)r (4.3) 



is assumed. With such W{i,j) the Hamiltonian ( |4.1| ) is a selfadjoint operator on the domain 
V{H) = © V^^ with = V{hi)®C. Note that, since H2 = C, the quantities W22(i - j) 

are complex numbers. The Wi2{i — j) are vectors in Hi, the W2i{i — j) are continuous linear 
forms on Hi, and the wii{i — j) are bounded operators in Hi. 

Let us denote by the Cartesian product = T x T x . . . x T of the N intervals 

N 

T = [— 7r,7r] and by p the points of T^, p = {pi,P2, ■■■,Pn), Pk ^ T. The assumption ( |4.2| ) 
implies that the operator-valued function 

n(p) = ^ WU) exp(i(p, j)), n{p) -.H^H, p G T^, (4.4) 
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where (p, j) = J2k=iPkjk is continuous and bounded on T^. Due to Eq. (|4.3| ), the values 

\UJ2l{P) ^22[P) J 

of this function represent selfadjoint operators in Ti for any p G T^, with 

^mn{p) ■ Hn^Hm] [t^ll(p)]* = t^ll(p), t^22(p) ^ R and [t^2l(p)]* = ^uip)- 

The quantity b{p) = u!i2{p) can be considered as a vector of Hi while co'2i(p) = ( •, b{p)) 
(cf. definition (|2.1|) of the Hamiltonian A). 

The blocks Rij,k;z), Rij,k;z) : ^ n^^\ j,k G Z^, of the resolvent R{z) = 

{H — z)^^ satisfy the equation 

iA-z)R{j,k;z)+ Yl W{j-f)Rif,k;z)=5,J. (4.5) 

After Fourier transformation in Q, 

{Fu){p) = ^^-^ ^^'■^exp(i(p,j)), p G T^, 



the system ( [4.5| ) takes the form 

(A - z)R{p, p'- z) + Vt{p) R{p, p'; z) = 6{p - p') I (4.6) 

where the quasi-momenta p, p' run through the set and R{p,p'; z) stands for the trans- 
formed resolvent R{z). Thus, the factorization ( ^.8| ) holds with 

G{p-z) = [A + n{p)-zr\ (4.7) 

First, let us consider the specific case where intercellular interactions W{i — j) in the 
Hamiltonian ( [4.1|) have the simple form 

W{t-j)= f ° r V (4.8) 

i. e., where the cells interact with each other only via the molecular channels. Obviously, in 
this case the factor G{p; z) is still given by p.9|). As compared to this expression the only 
difference is that now p G and 

M2 (p; z) = \2-z + UO22 (p) -I3{z). (4.9) 

Let 00^2^ = min c<j22 (p) and 00^^^ = max 0^22 (p) • Similarly to the analogous assumption 

in Sec. |^ we assume that the closed interval [A2 + ^^22^, X2 + ^^22^^] is totally embedded in the 
absolutely continuous spectrum adhi) of hi, and no thresholds of this spectrum belong to 
[A2 + uj^^, X2 + uj^^]. We also assume that this interval belongs to the holomorphy domain 
V of the function (3{z), and that for /i G [A2 + u;^''', A2 + uj^''''] the inequality (^J) holds. 
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Let us consider the time evolution of the system described by the Hamiltonian ( |4.1|) 
with the simple intercellular interactions ( [4 .81) subject to the above conditions. We start 
again from a pure molecular state ip = ipi ® ^2-, llv^nll ^ Qn-, = 1,2, with = and 
— llv^zll = 1- The probability Pmoi(v^,^) to find the system at a time t > in the 
molecular channel is given by the analogue of (|3.18|) . As in Sec. ^ one finds the quasi- 
momentum representation 



and relations (|3.21| ), (|3.22| ) and ( 3.24| ) are still valid with the only difference that instead 



of ( |3.23|) X2{p) is now of the form 

Hp) = + UJ22{P) , PGT^. 



According to Eqs. (|2.14| ) and ( p.l5|) the main terms of the roots Znuci(p) and Zmoi(p) of 



the function ( [4. 91) in the unphysical sheet read 

2:nucl(p) =zi- " , (4.10) 

M + UJ22{P) - Zi 

Zmo\{p) = A2 + UJ22{p) + ^ . (4.11) 

X2 + UJ22{P) - Zi 

The asymptotic equation (|3.24|) now implies 

Praoi{^,t)= [ dp\J{p,t)\'\Hp)\'= [ dpexp{-T^^\p)t}\Mp)\' + eit) (4.12) 



with 



T^^\p) = -2Imz^oi(p) = -21m. 



A2 + CJ22(P) - Zi' 

As in Sec. |^ the background term e{t) is small for any t > 0, e{t) = 0(||6|p) and \^{t)\ ^ 1. 

Thus, if the real part E^j^^ of the "nuclear" resonance zi belongs to the interval 
[A2 + uj^^, X2 + uj^^] then there are "molecular" states ip which decay via the "nuclear" 
channel with a rate as close as possible to the maximal value ( |2.16|) . In this case the com- 
ponents 932 (p) are localized in a close neighborhood of the manifold 

A2+^22(P) =4'^ (4.13) 

in the quasi- momentum space T^. In particular, if the initial state is prepared such 
that the component f2{p) is nonzero only for the quasi- momenta p lying in the domain 
IA2 + 0iJ22{p) — E^R^l < (5r^''/2 with some small 5 > 0, then one should only integrate over 
this domain in the integral ( [4.12|) . In such a case, under the condition (|2.8| ), a lower estimate 

for the decay rate F^'' is given, as in Sec. 0, by ^ ^ ^ — Thus, by varying 6 one can 



get a rate as close as possible to the maximum ( ^.16[ ) 



R 



Now, let us consider briefly the case of more general intercellular interactions W{i — j) 
where all the components uj^nnip), m,n = 1,2, of the matrix Q{p), p E , can be nontrivial. 
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In this case the component 6*22 (p; z) of the factor ([4.7|) is given again by Eqs. (|3.20|) and ( ^l9| 
However, the function P{z) in (^]9|) is to be replaced by the modified function 



(3{p; z) = {[h, + oonip) -z]-'[b + b{p)], [b + b{p)]) 

where b{p) = uji2{p)- We make the natural assumption that the direct intercellular interac- 
tions in nuclear channels Wii(i — j) are so weak that the term uJii{p) produces only a very 
small perturbation of the initial "nuclear" resonance Zi generated by the Hamiltonian hi. 
More precisely, we assume that the resonance zi (p) generated by the perturbed Hamiltonian 
hi + ujii{p) has the property 

- ^i| < riJ^ for any peT^ (4.14) 

and that no other resonances arise in the domain T>. Another natural assumption is that 
the strength of the interactions Wi2{i — j) and W2i{i — j), i 7^ J, between "nuclear" and 
"molecular" channels of different cells is much weaker than the one of a single cell. This 
is why we can assume ||/||&|| ^ 1 and the following features. The elements b{p) G Hi 
are such that for any p G T" the function z) admits an analytic continuation into the 
domain T) of the unphysical sheet and in "D a representation of the type ( ^.61 ) holds for 

P{p;z) = ^f^^ +r'ip;z), (4.15) 
zi{p) - z 

with an explicitly separated pole term a{p) / {zi (p) — z) and a holomorphic remainder 
P'^^^^p-jz). As in Sec. p B| we assume that \a{p)\ = Ca{p)\\b\\'^, and that for any p G 
the limiting procedure ( |2.7| ) is possible for Ca{p) while Ima(p) -C Rea(p). For the re- 
mainder I3^^^{p] z), p G T^, z eV, the same statements are assumed to be valid as for the 
function (3'^^^{z). Under these assumptions one can repeat almost literally the study of the 
probability P^oiiv^t) &s performed above in case of the interactions (|4.8| ). We find again 
that the asymptotics of Pjaoi{f,t) is given by ( [4.12| ) with 



T'Pip) = -2Im 



A2 + UJ22ip) - Zi{p) 



Let us denote by (Rea)max the maximal value of the function Rea(p) on the manifold 
( [4.13| ). It is obvious that if one prepares the initial pure "molecular" state ip such that its 



component (p2{p) is localized in a close neighborhood of the subset of the manifold ( [4. 131 ) 



where Rea(p) = (Rea)inax, then for the probability P^oi{y^,t) the main qualitative result 
remains practically the same as in case of the interactions ( [4.8| ). Namely, varying the support 
of the component (p2{p) in one can achieve a decay rate of the state ip as close as possible 
to the maximal value of the width r^^\p) in ( [4.13| ). The main term 4(Re a)niax/riJ'' of this 
value is again inversely proportional to the "nuclear" width F^"* . 
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